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Abstract 

The paper confirms the existence of a special configuration (among the infinitive 
number of a priori possible virial states) which a B stellar (Baryonic) component 
may assume inside a given D dark halo potential well. This satisfies the d’Alembert 
Principle of virtual works and its typical dimension works as a scale length (we call 
tidal radius) induced on the gravitational field of the bright component by the dark 
one. Its dynamic and thermodynamic properties are here analyzed in connection 
with the physical reason for the existence of the Fundamental Plane for ellipticals 
and, in general, for two-component virialized systems. The analysis is performed by 
using two-component models with two power-law density profiles and two homoge¬ 
neous cores. The outputs of this kind of models, at the special configuration, are 
summarized and compared with some observable scaling relations for pressure sup¬ 
ported ellipticals. The problem of extending the results to a general class of models 
with Zhao (1996) profiles, which are more suitable for an elliptical galaxy system, is 
also taken into account. The virial equilibrium stages of the two-component system 
have to occur after a previous violent relaxation phase. If the stellar B component is 
allowed to cool slowly its virial evolution consists of a sequence of contractions with 
enough time to rearrange the virial equilibrium after any step. The thermodynamic 
process during the dynamical evolution is so divided into a sequence of transforma¬ 
tions which are irreversible but occur between two quasi-equilibrium stages. Then, 
it is possible to assign: a mean temperature to the whole B component during this 
quasi-static sequence and the entropy variation between two consecutive virial steps. 
The analysis allows the conclusion that the induced scale length is a real confinement 
for the stellar system. This follows from the application of the 1° Thermodynamics 
Principle under the virial equilibrium constraint, by checking how larger configu¬ 
rations turn out to be forbidden, according to the 11° Thermodynamics Principle. 
The presence of this specific border on the space of the baryonic luminous com¬ 
ponent has to be regarded as the physical reason why a stellar galaxy belongs to 
the Fundamental Plane (FP) and why astrophysical objects, with a completely dif¬ 
ferent history and formation, but characterized by a tidal radius (as the globular 
clusters are) lie on the same FP. An other problem addressed is how this special 
configuration may be reached. This is strictly connected with the problem of the 
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end state of the collisionless stellar system after a violent relaxation phase. Even if 
degeneracy towards the initial conditions is present on the FP, the mechanic and 
thermodynamic properties of the special configuration suggest this state may be the 
best candidate for the beginning of the B component virial evolution, and also give 
a possible explanation for why an elliptical is not completely relaxed in respect to 
its dark halo. 

Key words: Celestial Mechanics, Stellar Dynamics; Galaxies: Clusters. 


1 Introduction 


As Ogorodnikov (1965) has highlighted, in order to find the most probable 
phase distribution function for a stellar system in a stationary state, the phase 
volume has to be truncated in both coordinate and velocity space. While in 
the velocity space the truncation arises spontaneously due to the existence of 
the velocity of escape, the introduction of a cut-off in the coordinate space 
appears, on one side, necessary in order to obtain a finite mass M and radius 
R , but, on the other, very problematic. 

A similar difficulty also appears on the thermodynamical side, for which an 
extensive literature exists ( from: Lynden-Bcll & Wood, 1968; Horowitz & 
Katz, 1978; White & Narayan, 1987, until, e.g., Bertin & Trenti, 2003, and 
references therein). By using the standard Boltzmann-Gibbs entropy: 

S = — J flnfd 3 xd 3 v (1) 


defined by the distribution function in the 6 — dimensional phase space , 
f(x,v) (hereafter DF ), and looking for what maximizes the entropy of the 
same stellar system, the conclusion is: the DF which plays this role in (1) is 
that of the isothermal sphere. But, the maximization of S, subject to fixed 
mass M and energy A, leads again to a DF that is incompatible with finite 
M and E ( see, e.g., Binney & Tremaine, 1987, Chapter 4; Merritt 1999, Lima 
Neto et al. 1999, Marquez et al. 2001, and references therein). 

It is beyond the scope of this paper to enter into the very complicate problem 
of looking for the suitable models for the collisionless stellar systems by an 
analysis in the phase space and a research of the DF which maximizes the (1), 
or to examine the thermodynamic properties of the family models which are 
able to explain the features of partially relaxed anisotropic stellar systems (see, 
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e.g., Stiavelli & Bertin, 1987, Bertin & Trenti, 2003, and references therein). 
Nevertheless, our limited contribution to the wide discussion existing in the 
literature will be to underline as in a stellar component, embedded in a second 
dark matter subsystem (as realistically thought, e.g., Ciotti, 1999, and refer¬ 
ences therein), a truncation is spontaneously introduced in coordinate space, 
due to the presence of a scale length induced from the dark halo, as long as 
virial equilibrium holds. That is the tidal radius which we discovered has to 
exist when two-component models are considered with two different power-law 
density distributions and two inner homogeneous cores (Secco, 2000; Secco, 
2001, hereafter LSI), under some constraints on the exponents. 

The consequence of the existence of a special configuration characterized by 
this tidal radius are analyzed here. 

Finally, we will gain more insight into the physical meaning of the special con¬ 
figuration considered by introducing the thermodynamic information quantity 
(Layzer, 1976). 

Even if some considerations which follow are more general and may also be 
extended to spirals, we will limit our considerations to the collisionless stellar 
systems, as the ellipticals are considered. 


2 Looking for a special virial configuration 


To introduce the problem in a general way, we start by considering the poten¬ 
tial well of a given spherical virialized dark matter (hereafter, DM) halo of 
mass Mr, and virial radius an, with a density radial profile as follows: 


p(r) 


_Po_ 

(■ r/r 0 p[l + (■ r/r 0 ) a ] s 


; S = (P- 7)/a 


( 2 ) 


where p Q and r 0 are its characteristic density and its scale radius, respectively. 
These kinds of profiles have already been introduced by Zhao (1996) and by 
Kravtsov et al. (1998) in order to generalize the universal profile proposed by 
Navarro, Frenk & White (hereafter, NFW) (Navarro et al.1996, Navarro et al. 
1997)which is obtained from eq.( 2 ) as soon as (a — 1 ; (3 — 3; 7 = 1 ; S — 2 ). 
Hereafter, we will name them Zhao profiles. 

The question which arises is the following: Does a special virial configura¬ 
tion exist among the infinitive number of a priori possible virial configurations 
which the luminous (Baryonic) component (B) may assume inside the given 
dark one (D)? 
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Fig. 1. Two-component model: D is the dark matter ( DM ) halo spheroidal compo¬ 
nent, B is the bright (Baryonic) inner one. Here the B component is non-homothetic 
to the D one instead of what occurs in this paper. Fig.l shows a B spheroid with 
an axis ratio smaller than that of the D one in order to show what is, in the general 
case, the dark matter fraction which, according to Newton’s first theorem, exerts 
dynamical effect on the embedded B subsystem. In this general case the dark mat¬ 
ter fraction is inside the surface S* which, in the homothetic case, coincides with 
the B contour (Raffaele, 2003). 

2.1 Tensor virial formalism 


In order to find the answer we need to use the tensor virial theorem extended 
to two components: D + B (Brosche et al. 1983; Cainimi et ah 1984; Caimmi 
& Secco, 1992). In a two-component system in which one (B) (Baryonic or 
Bright; in this context also: stellar) is completely embedded in the other (D) 
(of DM) and each of them is penetrated by the other, the following tensor 
virial stationary equations hold: 

2 (Tu)ij =-(V u )ij; (u = B,D ; i,j = x,y,z) (3) 


where ( T u )ij is the kinetic-energy tensor and (V u )ij is the Clausius’ virial tensor 
which splits into two terms: the self potential-energy tensor, (Q,,)^-, and the 
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tidal potential-energy tensor, (V uv )ij (u,v = B,D), due to the gravitational 
force which the v subsystem exerts on the u one. The eqs.(3) yield the following 
pair of tensor equations: 


2 (T B )ij - ~ {V BD )ij (4) 

2 {T D )ij = — (Vdb)p (5) 

Therefore, e.g., in the case of the inner B component, we have: 



where is the gravitational potential due to the mass distribution of the D 
component. 



Fig. 2. The energy trends of the 5-system as a function of size ratio of Baryonic to 
DM components, x = cib/cidj normalized at the factor (GM^F) / a B - The F-curves 
represent the Clausius’ virial energies, the 5-curves the corresponding total poten¬ 
tial energies, in the cases of Tab.l (Raffaele, 2003). 
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Even if we do not look for the distribution functions which correspond to 
the models considered in the next subsection, it may be useful to remem¬ 
ber the link between the tensor virial quantities and the phase space, given 
by considering the moment equations of the second order, in the coordinate 
space, of the 11° Jeans equations which, in turn, are the second order mo¬ 
ments, in the velocity space, of the Boltzmann equation (Binney & Tremaine, 
1987, Chap.4; Chandrasekhar, 1969, Chapt.2). It follows that in a system of 
collisionless particles with istantaneous x position and v velocity, character¬ 
ized by a distribution function f(x,v), with spatial density in phase space, 
v(x) = / f(x,v) dn, and mass density p{x) = m* u(x) (m* being the average 
stellar mass), the kinetic-energy tensors , TP, lip, are defined as: 


T • 

J-i] 

T' = - 

v 2 


\ J P v i v i = T ij + ; 

pTJiVj dx ; n ij = I pcrP dx ; 


( 7 ) 

( 8 ) 


that is by the mean of square velocities (ryuf), by the square of mean velocities 
( VjVj ) of streaming motions, and by the random square velocity components 
(of-), respectively. The general expression for the mean vpjj is as usual: 

vpJJ =-^f fv&j dx (9) 


In the case of one single component, the Clausius’ virial tensor, V lv matches 
the self-potential energy tensor, , that is: 


bp Up 



( 10 ) 


$p is the tensor potential (e.g., Chandrasekhar, 1969) defined as: 

— xJ 

with G the gravitational constant. It should be underlined that, in general, the 
presence of one external component causes the non-equality between the total 
potential energy tensors of the subsystems and their Clausius’ virial tensors, 
as occurs in the case of the single component (see, e.g., LSI). Moreover, we 
should remember that Clausius’ tensors are made only by forces and positions 
and then the mass fraction of the dark outer component which enters in the 
B Clausius’ tensor is only that which exerts dynamic effects on B, according 
to Newton’s first theorem. 


~g[ p(x ')——— 


x 
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2.2 Two-component models 


Then we need to model the two components (Fig.l). The models we consider 
here are the same as in LSI. That is they are built up of two homothetic 
similar strata spheroids with two power-laws and two different homogeneous 
cores in the central regions. It should be noted that they correspond, to a 
good extent, 1 to deal with two spheroids with smoothed profiles of this kind 
(spherical case 2 ): 


Pd = 

PoD 

1 1 / r \rl 1 

Cn — 

on 

( 11 ) 


1 + (— r 

y r oD > 


r oD 


Pb = 

PoB 

c B = 

a B 

( 12 ) 


1 + — r 

WoB J 


ToB 



C B and C B are the two concentrations of the two components. In the LSI 
models, r oB and r oD were the radii of the two different homogeneous cores, 
which typically assumed one tenth of the virial radii, a B and an, respectively. 
According to these, the concentrations in the smoothed profiles both become 
equal to ten. The density profiles of kind (11) and (12), have the advantage 
that they may be considered as a generalization of pseudo-isothermal profiles 
which, in turn, may be regarded as sub-cases of the more general Zhao profiles 
when: 7 = 0 ; (3 = a = b, d. 

But a realistic elliptical model has to be: e.g., a stellar component with a 
Hernquist (1990) (hereafter, Her) density profile and a dark halo with a cored 
or non-cored NFW profile (that means, according to eq.(2), respectively: a = 
1; (3 — 4; 7 = 1; S — 3 and a = 1; (3 — 3; 7 = 0; S — 3; for the cored NFW, 
a — 1; (3 — 3; 7 = 1; S — 2, for the non-cored NFW) (as in Marrno, 2003, where 
two homeoidally striated ellipsoids are considered). 

Then, the problem of transfering the outputs obtained with two cored power- 
law profiles (which also hold, to a good extent, for the models with smoothed 
profiles (11), (12)) to the more general class of models with Zhao profiles given 
by eq.( 2 ), is still open, even if some preliminary considerations will be made 
in the subsect. 4.5. 

Our aim is to try to explain some scaling relations for elliptical galaxies, which 
are essentially relationships among the exponents of the three quantities: 

1 Generally, the deviations are only of a few per cent and less than the 10% related 
to the Clausius’ minimum location and to its value, respectively. 

2 Even if the considerations which follow are more general, for the sake of simplicity, 
we will often limit ourselves to the spherical case without losing the validity of 
spheroidal case, which may be recovered simply by introducing a form factor, F, 
depending on the axis ratio (see, eq.(18)). 
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r e =the effective radius, 

I e = =niean effective surface brightness within r e 
(T 0 =the central projected velocity dispersion. 

The advantage of a simple cored power-law model is that it is able to extract, 
in a completely analytical way, some of these main correlations, highlight¬ 
ing the interplay of the parameters. Moreover, this preliminary analysis may 
also underline what are simply details in the model and what, on the con¬ 
trary, is stictly connected with the physical reason for the existence of a FP 
for two-component virialized systems. That allows us to open the road for a 
generalization of the present results. 


3 The special virial configuration 


As shown in the previous papers, a special configuration for the B compo¬ 
nent arises inside these kind of homothetic (see, LSI) models (for the non- 
homothetic, see, Secco, 2000), which describe its evolutive pattern, obtained 
by contraction inside a D one, under the assumption that D is at fixed size 
and shape. For the sake of simplicity, we will assume that the outer component 
is frozen, without considering this constraint to be too essential in order to 
determine the main features of the dynamic evolution we are dealing with ( 
see, LSI). The main reason for this assumption is indeed that the masses of 
the two components are not equal, the outer one being about ten times the 
inner one. As a consequence, tidal influences between the subsystems are not 
symmetric; the one acting from inner to the outer is actually weaker than the 
reverse (Caimmi & Secco 1992). Moreover, even if a contraction effect is in¬ 
duced by inner density distribution on the inner regions of outer halo, during 
the dynamical evolution, as already underlined by Barnes & White (1984) , 
that effect does not cause a dramatic modification of the outer mass distribu¬ 
tion, if there is a supernovae- driven outflow, according to N-body simulations 
(Lia et al. 2000). On the other side, models in which this constraint has been 
changed with some less stringent additional conditions (Caimmi, 1994), seem 
to come to the same conclusion. 

The special configuration appears because a maximum in the Clausius virial 
energy trend, as B contracts inside D, (and then a minimum in the kinetic 
energy) exists under the following constraints on the exponents: 

0<6<3;0<d<2^(6 + d)<5 (13) 


The total potential energy, (E pot ) B , on the contrary is always monothonic 



Table 1 

Physical parameters of the models considered in the figures 2, 5, 6. The common 
density profile parameters (eqs.(ll, 12)) are: Cb = Cd = 10, m = 8.5. Moreover 
xt = at/ap; for the definition of vqb', ^ud', Vy, see text. 


cases 

b 

d 

x t 

UDB 

VQD 

4 

1 ) 

0.0 

0.0 

0.389 

0.300 

0.300 

0.300 

2 ) 

0.0 

0.5 

0.346 

0.300 

0.312 

0.333 

3) 

0.5 

0.5 

0.361 

0.312 

0.312 

0.313 

4) 

1.5 

0.5 

0.419 

0.367 

0.312 

0.254 


(Fig.2). Indeed, by definition: 
( E, po t)B = ^b + Wbd 


(14) 


where Hr is the self-potential energy tensor trace and Wbd is the interaction- 
energy trace of the tensor defined as: 





(15) 


As already underlined, in general, ( Wbd )ij does not match V t] (see, Caimmi 
& Secco, 1992). 

By referring to the LSl-models with profiles of kind (11),(12), the special 
configuration (which appears under the constraints on the exponents of the 
density power-laws and by taking a frozen D subsystem, as seen in the same 
paper), corresponds to the dimension of the B component characterized by 
the following semimajor axis, called tidal radius: 

vnB 1 M b \ 3rd 
Vy (2 - d) M d ' aD 

where M and a are, respectively, the mass and the major semiaxis of the 
subsystem considered. The coefficient uqb enters the integral which gives the 
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self- potential energy tensor of the B component (see, LSI), so it is a function 
only of the b exponent; the coefficient u' v is defined in the following way: 


9 

u v — 2 [{ v b)m{vd)m 


-l 


C 


— (b+d) 
B 


(3-d)[5- {b + d)] 


(17) 


where (^ b ) m , ( z ' d)m are defined in LSI. u' v , which is a function of b and d 
and of the two concentrations Cb = Cb, enters in the definitions of the tidal 
tensor trace: 


t /■ / r,M bM d 

Vbd — ~v v G - F; 

or 


(18) 


where F is the form factor (see, Marmo & Secco, 2003, hereafter MS3) and 
Mb is the fraction of D matter exerting dynamical effect on B , according to 
Newton’s first theorem. To a good extent it is given by: 

Mb = M D (—?~ d (19) 


The same mass fraction normalized to M B becomes: 


m = 


M d 

M b 


a B )3 -d 
On 


( 20 ) 


Moreover we dehne the total mass inside the B-structure as: 
M^ ot = Mb + Mb = Mb( 1 + wi) 


( 21 ) 


These quantities also enter into the definition of the 5-Clausius’ virial tensor 
trace as follows: 




— VQB 


GM% , GM b Mb „ 
- r ~ v v - r 


Ob 


Ob 


( 22 ) 


or in the form normalized by the factor QF : 


vr / 2 -d a B M d 

Vb — - iAz m x ; x — —; m = —— 

x v ’ a D M b 


(23) 


Why this configuration of the B component inside the fixed dark matter po¬ 
tential well is so special, may be understood by considering, in the next section, 
its mechanical and thermodynamical properties. 


10 



4 Special configuration: mechanical and thermodynamical proper¬ 
ties 


4.1 Mechanical arguments 


In order to understand the full meaning of the tidal radius , which is defined 
by eq.(16), and to be able to consider the thermodynamical processes of the B 
component, we have to look at the physics related to Clausius’ virial energy, 
Vb- By definition, the Clausius’ virial tensor trace is given by: 


Qb — 
Vbd = 


r ^ 

/"5Z 

r =1 

3 


< 9 $ 


x, 


X r 


B 


dx r 

d$D 


dx r 


Vb — 12 b + Vbd 

(24) 

■ daffi = j Pb (r ~b ■ Ib) dxB 

(25) 

dx~h = J Pb{tb ■ fo) dxb; 

(26) 


where Jr and fo is the force per unit of bright mass due to self gravity and 
the dark matter gravity at the point fs, respectively. By definition, the work 
done by the self gravity forces, L s , in order to assemble the 5-elements from 
the infinity is given by the self- potential energy 125 and the work done by the 
tidal gravity forces, L t , in order to put the B component together with the 
D one from infinity through all the tidal distorsions (see, Caimmi & Secco, 
2004), is given by the tidal potential energy Vbd- 

Then a small variation SVr for a small displacement 5fs of all B points, has 
the following meaning: 

SV B = 5L S + 5L t (27) 


4-2 Small departures from virial equilibrium 


We will now consider what is the mathematical form of the Vb potential en¬ 
ergy variation as soon as the B inner system contracts or expands its initial 
volume S a of a small quantity A S Q . Following Chandrasekhar’s analysis (Chan¬ 
drasekhar, 1969, Chapter 2) the following holds: by definition, the Clausius 
virial is a global, integral parameter of an extrinsic attribute, that means of a 
quantity which is not intrinsic to the fluid element (like pressure or density) 
but something, we name F(x), which it assumes simply by virtue of its loca¬ 
tion such as the gravitational potential and its first derivative. The variation 
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of the integral: 


5 J p B F dx B 

So 


J p B F dx B - J p B F dx B 


So + ^So 


(28) 


when the istantaneously occupied volume by the fluid changes from S 0 to 
S 0 + A.S 0 by subjecting its boundary to the displacement £(x, t ) = x — x 0 , may 
be transformed into the integral over the unperturbed volume, that is: 


5 J p B F dx B = J p B AF dx B (29) 


where A F is the Lagrangian change in F consequent to the displacement £. 

The extension of this analysis to two-component systems has been performed 
(Caimmi & Secco, 2004) with the following result: 


SSl B — 


3 o 

PbJ 2 tr^-dxs (30) 


So 


3 3 


SVbd = - / PbJ2J2 d Xb 

1 t^\ Ai dx k dx r 


r =1 
3 


8x r 


\ ' c i b 

Pd 2 ^ C—— 


So 


Mo 

3 3 d 

PDT.T.z'k 


r c)t 

r = 1 


Mo 


k =1 r= 1 


: dx k 


(■ x r ?p^) dx D (31) 

OX r 


where the unperturbed volume of D-component is M 0 and £'(x, t ) is the 
amount of the perturbation in the point domain of the same component. 


Under the assumption of a frozen dark component, vanishes and the main 
result holds: 


5V b = 5L S + 5L t ~ 



3 

PbYI & 

r= 1 


dd> b 
dx r 


dx B 



(x r 


d$p x 
dx r 


dx B = 


{SV B ) ao 


(32) 


By definition of the tidal radius , which is the B dimension at the maximum 
of its Clausius virial energy, at frozen a B , the (SV B ) aD is stationary at a t 
(Fig.2)(see, LSI). Therefore, by moving of a virtual 3 displacement Sa B , from 

3 Here virtual means at frozen D. 
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as = a,t, we have: 


6L s + 6L t ~{5V B {at)) aD = 0 


(33) 


This means that the configuration at a t satisfys the d’Alembert Principle of 
virtual works (see also LSI). The physical reason is the following: if, e.g., B 
contracts, less dark matter enters inside the Sb surface, in the meanwhile 
the self gravity increases. The opposite occurs if B expands itself. Therefore, 
even if both forces are attractive, the works which correspond to them, for 
a virtual displacement, are of opposite signs (see, LSI). Therefore, the tidal 
radius configuration is an equilibrium configuration even if not stable because 
the total potential energy of B has not a minimum (Fig.2). 

The consequence of these mechanical arguments with the property of the tidal 
configuration to be able to distribute in about the equal parts the self- and the 
tidal- energies (see, next subsection and LSI), is to yield some outputs which 
are in good agreement with the corresponding observable scaling relations 
related to the elliptical galaxies, and in general to the existence of a FP for 
two-component virialized systems, as we have already highlighted in the past 
papers (LSI, MS3). 


4-3 Scaling relations at the special configuration 


Even if the physical explanation for the main features of FP for virialized 
structures and in particular for the dynamically hot ellipticals ( Bender et al., 
1992), is still an open question, the two main roads which are present in the 
literature (see, e.g., Renzini & Ciotti, 1993; Ciotti et al., 1996; Bertin et al., 
2002) may be summarized as: 

a) The tilt of the FP is due to the stellar population effect. 

b) The tilt is due to a non-homologous structural effect. The light profile is not 
an universal de Vaucoleurs profile but a Sersic profile (1968) which changes 
with the galaxy luminosity: the Sersic index n increases as the luminosity 
increases ( weak homology, Bertin et al., 2002). 

Neither a) nor b) seem to give the proper answer to the scaling relations 
problem; the population effect disappears in the K-band (Maraston, 1999), 
contrary to what is observed (Pahre et al., 1998). Therefore, a metallicity 
sequence of an old stellar population may at most fit the trend observed in the 
B-band (Gerhard et al., 2001). Moreover with the weak homology the observed 
tightness of the FP appears to be hard to explain (Bertin et al., 2002). 

According to theory performed in the papers: Secco (2000), LSI, MS3, we 
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tried to explain the tilt by assuming a strict homology which does not imply a 
constant ratio M B /L due to the presence of a dynamical effect caused by the 
scale length induced on the gravitational baryonic held from the dark matter 
halo distribution. 

We come back to the virial equations (3), in the trace form, for a B-component 
completely embedded in a dark halo. The FP we obtain (see, LSI): 


r e = c 2 c 1 V lfG 2 A = 2; B = -1 

_ 1 +fn 

(_t 2 — - 

F[n nB + mv' v \ 

seems, at a first sight, not substantially changed in respect to the one we get 
by using the virial equations for a single component system. 

But if we consider the trace of eq.(4), at the special configuration with the 
condition that the bright component is pressure supported (i.e., peculiar ki¬ 
netic energy T pec is dominant with respect to the rotational kinetic energy 
T r ot ) we obtain: 

\Mb < U 2 >~ f - Q B-VBD \ . ^ << ^ ( 36 ) 

^ \ Z / a B =at 


(34) 

(35) 


where < cr 2 > is the mean square velocity dispersion of the stars. By adding at 
a t the equipartition between the self- and tidal energy, the previous equation 
becomes: 


( \m b ^V \ 2 " d 

\v' v GM B M D F ) 


(37) 


which, instead of eq.(34) yields the following FP: 


r e 


r^j 


<Jo~ d aff* m - 2=3 M b 2 ~ d 


(38) 


That means: 


If - 


l ~M B 2 ~ d 


( 39 ) 
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4-4 Output vs. observables 


I) From the first equality we have: 


A 


2 

2 -d 


( 40 ) 


II) Moreover, clue to the two relationships which connect the FP coefficients, 
A,B with the observed tilt (see, e.g., Djorgovski & Santiago, 1993), that is: 


A = 

1+ " ( 41 ) 

I b = -t^ 

we immediately obtain the tilt of the FP: 


CX-t 


1 — d 
3-d 


(42) 


and the other coefficient: 


B = 


3-d 
2(2 - d) 


(43) 


only as functions of the dark matter distribution. First of all it shond be noted 
that the quantities A,B,oit which define the FP and its tilt are independent 
of m. That means the galaxies which belong to the FP may have a different 
fraction of baryonic matter in respect to the dark one, but their dark matter 
density profile must be the same. To probe the issue of the first identity we 
choose, e.g.: 

d = 0.5 =*► A = 1.33, B = -0.83 and a t = 0.20 

in good agreement with the observations in the B— band. 

In the range: d — 0 P- 1 =>• A = 1 -P 2 ; —B = 0.75 -Pi; cq = 0.33 P 0 in 
agreement with the data related to all bands (Marmo, 2003, Tab.1.1 and the 
references therein). 

This allows us to expect (see, LSI) that in other families of galaxies with 
dark matter halos of the same kind of elliptical galaxies, as, e.g., the spirals, a 
Fundamental Plane also has to exist which has the same A, B and a exponents, 
totally independent of a completely different luminous mass distribution. This 
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is in good agreement with the universal FP discovered by Burstein et al. 
(1997). 

Ill) It should be noted that the values of A, B , cq which define the FP as a 
whole, are not directly linked with the past cosmological conditions, but they 
are only a function of the dark matter distribution d. This is in agreement with 
what Djorgovski already noted (1992) on the basis of Gott & Rees (1975), 
Gunn (1987), Coles & Lucchin (1995), occurs for the scaling relations in a 
CDM scenario. Indeed, from a cosmological point of view the FP means the 
following relationship (Djorgovsky, 1992), which we also recover in LSI: 


2 n rec + 10 = A( 1 - n rec ) - B( 12a t + 4n rec + 8), 


n rec — effective spectral index of perturbations. 


IV) But in the projections of the FP on the coordinate planes the depen¬ 
dence on the cosmological spectral index appears via the parameter, 7', we 
introduced in LSI, as: 


1 _ 1 + 3ay ec (M) _ 5 + n rec 

y(M) ~ 3 “ 6 


(44) 


where, according to Gott & Rees (1975b) and Coles & Lucchin (1995, Chapts.14, 
15 ), a rec is the local slope of the CDM mass variance, cr 2 M , at recombination 
time t rec , given by: 


din (Jm [tree) 

arec = dlnM ’ 

If the total energy of the system is conserved during the transition from a 
maximum expansion phase to the virialization, with its re-distribution among 
the collisionless ingredients by a violent relaxation mechanism, the following 
dependences on m and Mb , for the three main quantities of FP, hold: 


r= ( y(3 ^ } ; R = l/ 7 ' 

(45) 

/ e ~ mV4; I = i = 2 1 ' ~ (3 - 

7 (3 - d) 

(46) 

sn,rS 1 (3 - d) - i 0 1 y - 1 

a 0 ~ m M b s = - y(3 _ d) iS = - y 

(47) 
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Therefore, the projected scaling relation, such as the Faber-Jackson relation 
(hereafter FJ), turns out to be: 


(3-d)—/ 


4 7 


L ~ rri (3-<*F(7'-i) ^'-no-d) 


(48) 


not only directly related to the dark matter distribution, via the exponent d, 
but also to the perturbation spectrum, via 7'. 

For a typical galaxy dark matter halo of Mp ~ 1O 11 M 0 , it turns out that 
7' = 2 (Gunn, 1987). If d = 0.5, we obtain: 

L ~ Mg 8 ; r e ~ M° B 5 => I e = L/2nr 2 e ~ Mg 0,2 

It should be noted that I e decreases as Mb increases as soon as: 

3 — d > 7'; 

that is d < 1. On this mass scale, the FJ becomes: 


L 


r^j 


m 


0.16 


<7 


3.2 

o 


If d — 1, then a t = 0 and as a consequence: 


L ~ 


a 


4 

O 


without dependence on m. 

V) Moreover on the Clausius’ minimum we have: 

(M “)< = M d 1 + 7| bjj) 

If b and d are universal, then the consequence is the proportionality of the two 
masses and then of the two ratios: 

~ L/M b (49) 


VI) From an observational point of view, we know (Cappellari et ah, 2004) 
that, for a sample of E and SO, either fast rotators or nonrotating giant ellip¬ 
ticals, the following tight correlation holds: 

M/L ~ a°' 8 (50) 
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From the previous scaling relations, we obtain: 


2a t V(M D ) 

M b /L ~ 


( 51 ) 


which on this dark matter scale yields an exponent of <r 0 exactly equal to 0.8 
without distinguishing between and L/Mb, according to eq.(49). 

That also has to be compared with Jprgensen’s (1999) value: 0.76 ± 0.08. It 
should to be noted, as in the relationship (51), that the dependence on the 
other factor m is completely negligible. Indeed it turns to be: m 0 ' 04 . 

VII) Another issue is: the ratio of the dark matter fraction over total mass 
inside the bright radius as- 

At a t , it becomes: 


( Mp \ 1 

v M «A 1 + ^(2 -d) 

That means that it only depends on the luminous and dark density profiles. 
If they are both universal for the galaxy family considered, this dark matter 
fraction has to be the same for all members. 

It should be noted that this result is independent of the total mass ratio, dark 
over bright, m. 

If the most probable value for d is around 0.5 and b ranges from 2 -4- 3 
(Jaffe 1983; Hernquist 1990), the most probable values for turn out to 

range from 0.57 A 0.80, which corresponds to log (= 0.37 A 0.69, with 
= ^-50 ^ = 2-5- The agreement with the Jprgensen’s histogram 

in Fig.5 (Jprgensen 1999), related to early type galaxies in the central part of 
the Coma cluster when the same IMF is assumed, appears to be very good. 

But one of the most important issues of the present theory appears to be 
related to the physical reason which does cause the tilt. 

In order to have the tilt we need to have the maximum of Clausius’ virial 
energy. This, in turn, requires to have the equipartition between the self- and 
the tidal- energy of the stellar component. By considering the derivative of Vb 
of eq.(23) in respect to x and according to eq.(19) and the constraints (13), 
we conclude that the dark matter mass has to increase steeper than ( clb/cld )• 
This means, in turn, that pr, has to decrease less than 1/r 2 at the border of 
the bright mass, in order that tidal energy may overcome the self-energy from 
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this border forwards. But if we enter deeper into the fine play of the exponents 
we are able to deduce a stricter constraint. 

Going back to the dependence of I e on the mass ratio m and M R , the message 
of the (46) is: I e depends on the cosmological history of the galaxies, on dark 
matter distribution and on the baryonic fraction which is inside. But the 
ratio L/Mb i that is the tilt , is totaly independent of the cosmic perturbation 
spectrum and of the mass ratio m. It turns out to depend only on the dark 
matter density profile. Indeed, if we look at 

L ~ I e r 2 m t+2r M t b +2R 


where the exponents satisfy the following relationships: 


i + 2r = 0 
I + 2R = 2/(3 - d) 

Therefore, it is clear how the ratio L/Mb loses its direct connection with the 
cosmology given by y'. Moreover in order to have the positive (=observed) tilt, 
we need: 


2/(3 — d) <1 =*► 0 < d < 1 

This condition is stricter than the necessary condition for the maximum: d < 2 
For d — 1 =>- the tilt disappears. 

For 2 > d > 1 ==>• the tilt appears but it is negative ( the opposite of that 
observed). 

Therefore, the slope of FP tells us a constraint on the density distribution 
of dark matter halo. To have a positive tilt we need the DM mass has to 
increase steeper than (ob/od ) 2 at the border of the bright mass. This means, 
in turn, that pr> has to decrease less than 1/r. If the contrary occurs (=>- p^ 
decreases faster than 1/r but less than 1/r 2 ) the tilt changes its sign. That 
immediately underlines a problem with a NFW density profile concerning the 
inner part of the halo. The debate is still open. From the theoretical side what 
appears relevant is the conclusion of a recent paper by Mticket & Hoeft (2003) 
in which the constraint on the exponent in the central dark halo region has to 
be: 0 < d < 0.5, obtained by using Jeans’ equations. From the observation side 
we underline a very strict limit for the exponent 7 (7 < 0.8), completely in 
disagreement with the majority of simulations (7 > 1), determined by fitting 
the rotation curves obtained with high resolution tecnique, for a sample of 
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spiral galaxies in which dark halos density profiles of Zhao kind have been 
used (Garrido, 2003). 


4-5 On the limits of the m.odel 


The problem of generalizing the results obtained in the case of a two-component 
model with two power-law density profiles and with two homogeneous cores 
is still open. Some considerations may help in this future work. 

A) The existence of Newton’s first theorem (see, Fig.l), which allows us to 
ignore the mass distribution of the dark halo outside the B volume as re¬ 
gards the dynamical effects on the stellar subsystem and then for the trend of 
Clausius’ virial trace tensor Vb, which is the key of the whole theory. 

B) The conclusion of the previous subsection that is the presence of the tilt 
and its sign depends on the gradient of dark matter distribution at the border 
of bright mass confinement. 

Both together seem to converge on the main role which the central region of 
the dark halos plays in the present theory. That may occur in a satisfactory 
way only if the scale radius of the dark halo, r 0 o is not too small in respect 
to the dimension which contains the most of the stellar mass. In other words, 
we will expect that the dark halo concentrations have to be not too high. 

We show in Fig.3 and Fig.4 (Marrno, 2003) the same as in Fig.2 in the case of 
more general models which belong to Zhao models. That in order to prove as 
the appearence of Clausius’ minimum, for suitable values of the concentrations, 
seems a common feature of the these most general and most realistic triaxial 
models where the density profiles are given by eq.(2) for both the system 
components (e.g., NFW+Her, Fig.3; cored NFW+ Her, Fig.4). But an other 
problem arises, with this kind of models. That is to recover, analytically, the 
interplay of the exponents which appear in FP-quantities. 


4-6 Thermodynamic arguments 


We will now consider the B component thermodynamics inside a two-component 
system, to look for what characterizes its special virial configuration from the 
thermodynamical point of view. We begin to analyze the double system when 
it arrives at virialized stages after a phase of violent relaxation ( Binney & 
Tremaine, 1987). From this time onwards, it may be assumed that the B 
component begins its virial evolution which consists of a sequence of slow con¬ 
tractions witht enough time to rearrange the virial equilibrium after any step 
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0.25 0.5 0.75 1 1.25 1.5 1.75 2 


x 

Fig. 3. The trends of the Clausius’ virial energies in the normalized form, vs. x 
defined as in Fig.2, (dashed curves), in the case of two-component systems built up 
of a inner stellar component with a Hernquist (1990) density profile and an outer 
dark halo with NFW profile (eq.(2)). The concentrations are: Cb = 4 and, from 
top to bottom: Cb = 4, 7, 10. The mass ratio is m = 12. The curve corresponding 
to :Cb = Cb = 4, exibits surely a minimum for the absolute value of Clausius’ virial 
energy. For comparison, the self potential energy of the same B component when 
single, is shown with a continuous track (Marmo, 2003). 

of the sequence. In this way the thermodynamic process of contraction may 
be divided into a sequence of transformations which are irreversible but occur 
between quasi-equilibrium stages with the typical character of external ther¬ 
mal irreversibility (Zemansky, 1968). Therefore, it is possible to assign a mean 
temperature T s to the whole component during this quasi-static sequence of 
its dynamic evolution in the following way. 
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Fig. 4. The trends of the Clausius’ virial energies in the normalized form, vs. 
x defined as in Fig.2, (dashed curves), in the case of two-component systems 
built up of a inner stellar component with a Hernquist (1990) density profile 
(a = 1; fi = 4; 7 = 1; <5 = 3) and an outer dark halo with a cored NFW profile 
(eq.( 2 );a = 1 ; f3 = 3 ;7 = 0 ; <5 = 3). The concentrations are: C B = 4 and, from top 
to bottom: Cb = 4, 7, 10. The mass ratio is m = 12. For comparison, the self 
potential energy of the same B component when single, is shown with a continuous 
track (Marrno, 2003). 


By assuming that the B stellar system has an isotropic velocity distribution, 
we may define its mean temperature as: 


m* < cr 2 > 
= k 


(52) 


where m* is the mean mass of the stars and k is the Boltzmann constant (see, 
e.g., Lima Neto et al. 1999; the a -Ansatz, in Bertin & Trenti, 2003). 

Moreover, if we limit ourselves to the macroscopic pressure supported elliptical 
systems, which are relevant in order to define the FP we are dealing with, it 
follows that Ts is related to the dominant peculiar kinetic energy, \M b < 
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Fig. 5. Trends of the entropy function F(x) (eq.(66)) for the B-system, in arbitrary 
units (see text), as function of x (see, Fig.2). The cases are of Tab.l. The maxima 
occur at the tidal radii (Raffaele, 2003). 

o 2 Tbi in the following way: 


T 


s — 


2 T b 
N k 


(53) 


N being the star number of the B component. T B is, in turn, connected with 
the central projected velocity dispersion <y Q by the usual factor k v which links 
the kinematic galactic structure with a 0 as follows (e.g., LSI): 

o 2 0 = k v < a 2 > (54) 


In the 1° Thermodynamic Principle equation for the B system what has to 
appear is the work done against its pressure by both the forces on the system: 
the self gravity and the gravity which the dark matter distribution exerts on 
it. The corresponding potential energy variation, for a small contraction, will 
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be A Vb in such a way as to have: 


A E t = A Q- A V b 


( 55 ) 


where Et is the internal energy of B and A Q is the small amount of heat 
the structure is able to exchange with the surrounding medium, in which 
the two-component system lies,e.g., by cooling processes. Indeed, we have to 
take into account that the quasi-static evolution of the inner component ( 
the D one is considered non-dissipative) needs to involve some dissipative 
processes, as Prigogine (1988) has well pointed out, on the general grounds. 
In a gas dominated structure, dissipation may easily occur, e.g., by gas clouds 
collisions, but also in a collisionlcss stellar component such as an elliptical. 
Indeed, several evolutionary processes may potentially occur (see, e.g., Bertin 
& Trenti, 2003), some of them with a dissipative character. 

If enough time is available in order to restore virial equilibrium, the variation of 
the virial quantities of the same component during this quasi-static transition, 
turns out to be: 

A T b = -AVb/2 (56) 


Combining the two equations (55, 56) and considering that the internal energy 
may now be identified with the total macroscopic kinetic energy of the stars, 
Tb, according to eq.(53), we obtain the two equations: 


A Q = AVb/2 
AT 5 ~ — AVb/2 


(57) 

(58) 


The two requests that the system settles in virial equilibrium and obeys the 
1° Principle yields the equipartition of the virial energy variation which now 
becomes the variation of the Clausius’ virial energy: that is one half of it has 
to be exchanged with rest of the univese, the other half has to contribute to 
change the temperature of the system. The result is formally the classical one 
found by Chandrasekhar (1939) and Scwarzschild (1958) for a single gaseous 
stellar component. The big difference, for the B component inside the D one, 
is that the potential energy Vb is now a non-monotonic function of cib (Fig.2). 

The variation of the entropy of the B system, S b , during the transformation 
between two virial states due to a small contraction A as, which has the typical 
character of external thermal irreversibility, is evaluated as: 

AC! 1 AC'I o Ar , A (Vb/2) Szn\ 

AS B = =-A Q = 2 Nk——— (59) 

is —Vb/^ 
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If A V B is negative clue to the energy lost by radiation, the consequence of the 
radiation flow is an entropy amount AS r of the thermal radiation bath in the 
surroundings, which has the mean temperature T r , equal to: 


AS r 


1 (-AVb) 

2 T r 


(60) 


AS r is of opposite sign of AS^ and greater than it, in absolute value, in a 
way that, according to the 11° Thermodynamics Principle, the whole universe 
increases its entropy S u of: 

AS„=i(-Ar B )A-=i-) ( 61 ) 

Z ± r 1 s 


That holds because the surrounding thermal bath, in which the radiation flow 
will at the end thermalize, has a mean temperature lower than that of the 
stellar structure. 

Due to the non-monotonic character of the V B trend, these very important 
consequences follow: 

A) as soon as the Clausius’ virial energy is stationary S B must also be so. That 
occurs at at because the corresponding configuration satisfies the d’Alembert 
Principle of virtual works (see, eq.(33)). Therefore, one expects to obtain a 
maximum or a minimum for the entropy of the B component when the stellar 
system is on its tidal radius configuration. It should be underlined that this 
result is independent of the special class of two-component models we are 
dealing with but it depends only on the meaning of Clausius’ virial energy 
and on presence or not of the Clausius’ virial maximum during the quasi-static 
contraction sequence the chosen models are able to describe the B subsystem 
evolution. Indeed, the eq.(59) derives from the physical meaning of Clausius’ 
energy and from the first Thermodynamical Principle together with the virial 
constraint. 

B) An other consequence of eq.(59) is: all the configurations of the stellar 
system which correspond to a dimension greater than the tidal radius at are 
forbidden. Indeeed let us consider any configuration on the right side of the x t 
in Fig. 2. Starting from this one we take into account the thermodynamical 
transformation which follows to a small contraction A a# < 0. According to 
the eqs.(53, 57, 58, 59), the consequences are: 

A a B < 0; A T s <0; AQ > 0; AS B > 0 (62) 


The stellar structure would have to increase its entropy by taking energy from 
the radiation bath in which it is embedded and which has a lower temperature 
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than Ts■ It is manifest that this termodynamical process contradicts the sec¬ 
ond Principle. Starting from the same configuration, we now consider a small 
expansion. We would obtain: 

A a B > 0; A T s >0; A Q < 0; AS B < 0 (63) 

From the thermodynamical second Principle that would be possible but the 
contradiction arises from the conservation of the total energy which is equal 
to the half ot the total potential energy, \E pt , of the two-component system 
(see, MS3, Fig.5). Indeed a small expansion would request an increase of the 
total mechanical energy of the whole system. 

The main result is that all the virial configurations greater than a t are forbid¬ 
den. 

C) On the contrary, all the configurations characterized by a semimajor axis 
a B < (k correspond to obtain, for a small contraction, the transformations: 

A a B < 0; A T s >0; A Q < 0; AS B < 0 (64) 

which are allowed by the second Principle. A small expansion would on the 
contrary yield to: 

A a B > 0; A T s <0; A Q> 0; AS B > 0 (65) 

which are forbidden both from the energy conservation Principle and from the 
second Thermodynamical Principle. The last pair of relations are those we may 
obtain by using the eqs.(57, 58, 59) in the case of one single B component. That 
is clearly because: for a B < a t the regime is characterized by the overcoming 
of self-gravity on the tidal-gravity. 

The main conclusion is: the scale length induced on the stellar component by 
the dark matter halo works as a real border of this subsystem in the same 
way as the Hoerner’s (1958) tidal radius, induced by the galaxy, acts as a 
confinement for the stellar of a globular cluster (see, Appendix A). Therefore, 
it appears rasonable that ellipticals and GCs belong to the same FP and that 
King’s (1966) models, which need tidal cut-off, could be shared by both kinds 
of objects (Djorgovski, 1995; Burstein et al. 1997; Secco, 2003). 

Moreover, the possible way in which the special tidal radius configuration may 
be reached during the violent relaxation phase of the system, gains a deep 
meaning. Indeed, this configuration is the widest one, the baryonic matter 
may have inside the dark potential well, and corresponds to the minimum of 
the macroscopic random velocity pressure the stellar system has to gain in 
order to virialize itself. 
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5 Entropy trend and thermodynamic information 


By integration of eq.(59) (regarding the small variations as infinitesimals) we 
may obtain the trend of the entropy, normalized to the factor 2 Nk, S(x), 
along the quasi-static virial sequence of the B component, as follows (see also, 
MS3): 

S(x) — S(l) = In | = Fix)] x = ob/od (66) 

V B (x) 

where x — 1 is obtained when both the two subsystems coincide. The absolute 
value of Vb(x) exibits a minimum at x = Xt along the virial sequence, therefore 
the function F(x) has its maximum at the special configuration of B. The 
corresponding entropy variation, in physical units, along the sequence is simply 
2NkF(x). 

In the case of similar heterogeneous models considered, for the four cases of 
Tab.l, we obtain the trends shown in Fig.5. 

The last thermodynamical quantity to be considered will be the thermody¬ 
namic information, which, according to Layzer (1976), is: 

I = Smax - s (67) 

where S ma3 , means the maximum value the entropy of the system may have 
as soon as the constraints on it, which fix the actual value of its entropy to 
S, are relaxed. Therefore, in order to increase its information a system has 
to decrease its entropy more and more in respect to that of the universe (the 
maximum available and given essentially by the entropy of CBR ; see, e.g.: 
Coles & Lucchin, 1995; Secco 1999). 

For a luminous component such as B, which is embedded in an other D, a 
contraction decreases the entropy function F(x) only if its dimension is smaller 
or equal to a t . If we now transform the information given by the eq.(67), into 
the same units of F(x ), we obtain: 

Imin = Smaz - S(l) - F(x t ) (68) 

This means that from tidal radius downwards the component B may gain 
information by contraction. We plot in Fig.6 the trends of the information 
/ in the same units of F(x) by assuming S max = S(ay) instead of the real 
maximum entropy value corresponding to the CBR. The eq.(67) becomes: 

7 = S(x t ) - S(l) - F{x) (69) 
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Fig. 6. Thermodynamical Information trends for the B- system, in arbitrary units 
(see text), as a function of x (see Fig. 2) for the cases of Tab.l. The maximum 
normalized value S max is here assumed to be S max = S (xt) instead of the real 
maximum entropy value corresponding to the cosmic background radiation (see, 
eq.(69)) (Raffaele, 2003). 


The tidal radius configuration appears again as the best candidate for the be¬ 
ginning of the virial stage because it corresponds to the minimum of thermo¬ 
dynamical information along the whole virial sequence. The stellar subsystem 
has the real possibility to evolve onwards, with the dissipative processes indi¬ 
cated, e.g., by Bertin & Trenti (2003), by becoming more structured than it 
does at the end of the relaxation phase, that is at the special configuration 
given by x t . The same arguments may be transferred to the cut-off spirals 
which, in this frame, had to begin a significant structure evolution from this 
stage forwards. Even if this extension has to be made wider and deeper in the 
future, it may become the ground for the interpretation of the cut-off radius 
observed in many edge-on spiral galaxies (see, e.g., van der Kruit, 1979; Pohlen 
et al. 2000a, 2000b and Kregel et ah, 2002). A preliminary analysis has been 
done in Guarise et al. (2001) and in Secco & Guarise (2001). 
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6 Concluding remarks on the special configuration 


The consequence of the existence of a special configuration characterized by 
the tidal radius have been analyzed. We stress particularly that: 

a) at this special configuration the two-component models considered, built 
up of two homothetic similar strata spheroids with two power-laws and two 
different homogeneous cores, yield outputs which are relevant for some observ¬ 
able scaling relations of pressure supported ellipticals (subsect. 4.3 and 4.4). 
Moreover they could explain the physical reason for the existence of a Funda¬ 
mental Plane for two-component virialized systems (e.g.,Dressier et al.1987; 
Djorgovski & Devis, 1987; Bender et al. 1992; Djorgovski & Santiago, 1993; 
Burstcin et al. 1997; Bertin, Ciotti & Del Principe, 2002; Borricllo, Salucci & 
Danese, 2003, and references therein). 

b) The models belong, once their profiles are smoothed, to the class with pro¬ 
files of general pseudo-isothermal kind, in turn, a sub-case of the more general 
class with profiles introduced by Zhao (1996). Indeed the density distributions 
of a realistic two-component model for an elliptical belong to this last class: 
both the dark matter halo profile, as proposed by Navarro, Frenk & White 
(Navarro et al. 1996, Navarro et al. 1997) (see section 2) and the Hernquist 
(1990) profile for the stellar component. 

Even if the analytical properties we consider here are those we deduce by using 
the non-smoothed power-law profiles with two homogeneous cores, they may 
be transfered, to a good extent, to the general pseudo-isothermal case. The 
other problem, to recover analytical, similar results also in the more realistic 
model with, e.g., NFW+Her profiles, is still open and has been adressed in 
the subsect. 4.5. 

c) The thermodynamical relevance of this induced scale length is high in con¬ 
nection with the Fundamental Plane of pressure supported ellipticals. Indeed, 
from the point of view of the two-component galaxy thermodynamics, deduced 
by the 1° Thermodynamic Principle under the virial equilibrium constraint, 
this dimension works like a wall of a vessel, because a larger configuration turns 
out to be forbidden by the 11° Principle of the Thermodynamics. This cut-off 
on the luminous component space provides the gravitational field, which is in¬ 
trinsically without any scale length, of a specific border, as it appears for other 
stellar systems such as the globular clusters. As we have already shown (LSI, 
Secco, 2003; MS3) this truncation, which King (1966) has also introduced ad 
hoc in his primordial models for ellipticals, seems to be the common feature 
of all the astrophysical structures, from the globular clusters to the galaxy 
clusters, belonging to the cosmic metaplane, (Djorgovski, 1995; Burstein et 
al. 1997). This is in the former induced by the Galaxy potential well, in the 
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others by the dark matter halos. Its existence might be able, in principle, 
to explain some of the fundamental scaling relations of these structures as 
we already have proved for the ellipticals (LSI) and, tentatively, for the GCs 
(Secco, 2003). How this special configuration may be reached, during the stel¬ 
lar system evolution, is strictly connected with the problem of the end state 
of the collisionlcss stellar system after a violent relaxation phase and then to 
the problem of the constraints under which this phase occurs (Merritt, 1999). 
What we have already proposed in LSI, is that the end state has to corre¬ 
spond to the minimum of the macroscopic pressure that the stellar system 
needs in order to virialize. This is given by the maximum of Clausius’ virial 
energy configuration which, we will prove, also has the property of maximiz¬ 
ing the entropy of the stellar component located in a two-component virialized 
system. It corresponds to the wider configuration which does not break the 
second Thermodynamic Principle together with energy conservation. This en¬ 
forces the idea of looking at this special configuration as the best candidate 
for the end of a violent relaxation process of a stellar system when it occurs 
inside a dark massive halo component. Indeed, the presence of this widest al¬ 
lowed configuration for the baryonic matter, which has the least requests for 
sustaining the structure in virial equilibrium, would justify the fact that the 
ellipticals are not completely relaxed systems in respect to the collisionless 
dark halo. As White & Narayan (1987) have pointed out, by studying single 
power-law stellar structures, the ellipticals seem indeed to have stopped their 
violent relaxation process before its end, unlike the collisionless dark matter 
structures. Really, the NFW profile is given by eq. (2) with the exponent (3 = 3 
instead of the Hernquist profile, which is in agrement with the de Vaucouleurs 
light profile, and requires (3 = 4. 

Moreover, it should be underlined that, due to the physical meaning of Clau¬ 
sius’ virial energy, the maximization of entropy inside the virial evolution 
sequence is guaranteed as soon as the model exibits the minimum of the ab¬ 
solute value of Clausius’ virial, apart from the specific model we are dealing 
with (see, subsect. 4.6). 


7 Conclusions 


The Mechanics and Thermodynamics of a stellar virialized system, embedded 
in a dark matter halo, have been considered. The models (the same used in 
LSI) are not derived by looking for the DF which maximizes the standard 
Boltzmann- Gibbs entropy. They consist of two heterogeneous spheroids with 
two power-law profiles and with two homogeneous cores. The concentrations 
are both equal 10. Linder some restrictions, they may be considered as belong¬ 
ing to a sub-class of the general class of models endowed with Zhao density 
profiles. The most relevant aspects are: 
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(«) 

- Under the usual constraints on the exponents of the density profiles and the 
assumption of a frozen dark matter halo, a configuration exists, characterized 
by having a dimension called tidal radius (at), which enjoys special properties 
both from the mechanic and the thermodynamic point of view. 

The mechanic properties at (^-configuration are: 

i) -the absolute value of the Clausius’ virial energy reaches its minimum. 

ii) -The d’Alembert Principle of virtual works is satisfied, therefore it is a con¬ 
figuration of equilibrium. Nevertheless it is not of stable equilibrium because it 
does not minimize the total potential energy (E pot ) B of the stellar subsystem 
(see, Fig.2 and MS3). 

This special size translates to FP-like relation by yielding outputs very rele¬ 
vant for some observable scaling relations connected with pressure supported 

ellipticals. 

To fit the observed FP would require the dark halo density distribution to 
be close to a generalized isothermal profile with an exponent d around the 
value 0.5 at the border of Baryonic mass confinement. A stringent upper limit 
appears to be d = 1 for which the gradient of dark matter distribution would 
produce a tilt exponent: a t = 0. That means: the observed tilt disappears. 

The thermodynamical properties at (^-configuration are the following: 

A) -The last mechanical property has the consequence, from the thermody¬ 
namical point of view, that the work done from the external forces, i.e. the 
self gravity and the tidal gravity, against the macroscopic pressure of the B 
stellar subsystem, is zero at at for any small contraction or expansion starting 
from it. This also means that the entropy of B is stationary at a t , if one consid¬ 
ers the first Thermodynamics Principle under the virial equilibrium constraint 
for the thermodynamical irreversible transformations occurring between two 
consecutive quasi-equilibrium stages in which the evolutive sequence may be 
divided. At every quasi-equilibrium configuration it is possible to associate a 
mean temperature and at two consecutive stages, an entropy variation. 

B) -By integration of the entropy variation, it appears that entropy reaches, 
inside the evolutive virial sequence, its maximum at a*. 

Therefore, the first conclusion is that: as soon as a class of models is able 
to exibit a minimun of the virial energy along an evolutive sequence of the 
two-component structure they describe, the model of this class which corre¬ 
sponds to the Clausius’ minimum (in absolute value) is also the model which 
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maximizes the entropy inside the virial sequence. This is verifyed in spite of 
the class of models we are dealing with and it appears new in respect to the 
results, for single structures, which appear in the literature. Moreover, this 
also allows us to avoid a pragmatic approach in which the thermodynamic 
request is separated from the mechanical one as in Lima Neto et al. (1999). 
Indeed the configuration at at enjoys both the best properties, mechanic or 
t her mo dynamic. 


iP) 

- The second relevant aspect is that the mechanic and thermodynamic proper¬ 
ties together reveal the confinement capability of the scale length induced on 
the stellar subsystem from the dark matter halo, That turns out to have the 
same role of the tidal radius induced by the Galaxy on a globular cluster, as 
Hoerner (1958) found, and represents a generalization of his result (see, Ap¬ 
pendix A). This seems to be the key for understanding the physical reason for 
the FP structure and why objects with a completely different history of forma¬ 
tion and evolution share the same FP features. Indeed, the second Principle 
of Thermodynamics and the energy conservation transform this tidal radius in 
a true cut-off, i.e. the stellar system, as a whole, cannot have a configuration 
wider than that which corresponds to the special one. That may be regarded 
as connected with the request, coming from many open problems, as referred 
in the introduction, of a tidal cut-off on the coordinate space for the DF. 


(7) 

-The third relevant aspect is related to the possible way in which this special 
configuration may be reached. Since the beginning (LSI) we noted that Clau¬ 
sius virial maximum corresponds to the minimum of the macroscopic pressure 
a subsystem needs for virialize during the relaxation phase and the consequent 
conversion of radial ordered velocity into a dispersion velocity field (see, e.g., 
Huss et al. 1999). This now becomes more than an Ansatz. Indeed, it ap¬ 
pears that the widest allowed configuration for the stellar system corresponds 
to this minimum. This may explain how and why this special configuration 
will be reached by solving also the problem adressed by White & Narayan 
(1987). A possible physical reason is the existence of this (^-configuration for 
the baryonic component, needed because it is less massive than the dark halo 
and therefore more affected by the other component tidal effect than the halo 
does. 

On the other hands the problem of knowing exactly the end state of a col¬ 
lisionless stellar system after a violent relaxation phase, would require the 
knowledge of the initial conditions under which this mechanism occurs (Mer¬ 
ritt, 1999). But as Djorgovski (1992) has pointed out and we have shown in 
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LSI 4 a degeneracy exists on the FP towards the cosmological initial condi¬ 
tions, then the possible end state of a stellar system has probably to be only 
indirectly deduced from the features of FP. The capability of Clausius virial 
maximum configuration to explain some of the main features of the FP, may 
highlight what is, otherwise, covered by intrinsic FP degeneracy. 


(^) 

-The fourth is: taking into account the thermodynamical information during 
the evolution sequence considered here for the B stellar system, it arises that 
the special configuration at a t corresponds to the minimum of this informa¬ 
tion. That means an elliptical will gain small information due to its small 
dissipation capability and therefore will stay at a low level of structuration. 
On the contrary, a gas dominated galaxy will become more structured and dif¬ 
ferentiated by this configuration forwards due to its dissipation engine given, 
e.g., by cloud collisions. 
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A Appendix 


On the Hoerner’s tidal radius 


How is the relationship between the tidal radius an induced by the Galaxy 
on the satellites, as the Globular Clusters (hereafter GCs), which Hoerner 
(1958) discovered and the special size which maximizes the Clausius’ virial 
energy of the Baryonic component? The roles of confinement are the same 
(see, e.g., (/^-statements of Conclusions) and the mathematical structures are 


4 The degeneracy disappears for the ratio of the dark over baryonic mass , rri, (see, 
LSI). 
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very similar. The Hoerner’s result is indeed: 



where, M c , is the GC mass and M G is Galaxy mass contained in the a G ra¬ 
dius. The eq.(A.l) has to be compared with the general result obtained in the 
contrentic case given by eq.(16). a# is obtained by balancing of all the accel¬ 
erations on the Cluster point P in Fig.lA, when the variation of centrifugal 
force, due to a rigid Galaxy rotation, at the point P is taken as negligible in 
respect to that of the Cluster center O (see, e.g., Brosche et ah 1999). It should 
be underlined that the eq.(A.l) is derived in the Hoerner’s point mass model. 
On the contrary, the eq.(16) is deduced by the non monotonic character of the 
Clausius’ energy dne to a mass extension of a component submitted both to 
its self gravity and embedded in the gravity of an other one. 

The balance to consider is not of forces but between self- and tidal- potential 
energy with the same sign but with different trend. The question is: in the 
concentric case only comparison between energy may be possible (the forces 
are concordant), but in the off-center case both the definitions are possible? 
That means: in this last case also the Clausius ’virial of the Cluster shows a 
non-monotonic character? 

By analogy with the concentric case, we already guessed (Secco, 2003) that, 
for some mass distributions of Cluster and Galaxy, the answers were positive. 
Detailed computations of the self- and the tidal-energy tensors have been 
already performed in a previous paper (Caimrni & Secco, 2003) where the 
general second order theory in a c /R 0 ( R a is the mean orbit radius of the CG) 
for the two off-center component systems is developed. For heuristic sake, 
we take now into account the simple case in which both the two off-center 
components, Cluster and Galaxy, are considered as two homogeneous spherical 
mass distributions (Cainnni & Secco, 2005). 

Generally speaking, the global tidal-potential energy trace due to the dynam¬ 
ical effect of the Galaxy-matter distribution on the cluster, V cG , may be split 
into: 


Kg = Vf G + K°g 


(A.2) 


where Vf G is the potential energy of a mass point placed at the cluster barycen- 
tre, with same mass as the cluster, due to the Galaxy mass fraction Ad G (R 0 ) 
which is located inside the mean cluster distance R 0 , and V' cG represents an 
additional term which is due to the cluster mass distribution. The same holds 
for the global kinetic energy, T c , which, according to Koenig’s theorem, may 
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be split into two terms: 


T c = T' + T° 


(A.3) 


one term, Th is the intrinsic kinetic energy of the cluster, the other one, T°, is 
the kinetic energy of the cluster barycentre. Because of the time mean motion 
of the cluster barycentre occurs in stationary virial equilibrium, the following 
relation holds: 

2 T° = -V° G (A.4) 


as soon as the two quantities are averaged over the period of cluster orbit 
inside the gravitational Galaxy held. Therefore, the virial equilibrium of the 
cluster mass distribution gives: 

2 Z = -K = -a, - Vg (A.5) 


where V' c is the Clausius’ virial energy of the cluster due to its mass distribution 
and fl c is its self-potential energy. 

From the second order theory in a c /R 0 , the following result may be deduced: 
Vg = -ic—(^) 3 (A.6) 

O Qj C IXq 


Then the Clausius’ virial related to the mass distribution of the globular cluster 
becomes: 


3 GMl _ \ c M c M g {R 0 ) ,a c ^> 
5 a c 3 a c ^ R a ' 


(A.7) 


The non-monotonic character of V' c is manifest and then it follows that also 
in the case of two off-center component system, the Clausius’ virial exibits a 
maximum given by: 


o>Gt — 


/9\i/3/l M c 

vy V2 Mg 


ac 


(A.8) 


which turns out to be only about a factor 1.2 greater than the a H radius. 
The corresponding thermodynamical arguments in the case of two off-center 
component system are still an open question. 
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Fig. A.l. Two off-center spherical component model. The satellite (e.g., a Globular 
Cluster) is embedded in the gravity field of an other subsystem (e.g., the Galaxy). 
Its mass M c is distributed inside a c and its barycenter O moves in a mean circular 
orbit (dashed circle) of radius R 0 around the center O' of the other component with 
mass Mq inside the radius ac- P is the point on which the accelerations hold the 
balance when a c becomes the tidal radius of Hoerner. 

B Appendix 


Glossary of Simbols (in order of comparison) 
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self potential-energy tensor 
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self potential-energy tensor 
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Section 5 
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IX — S thermodynamical information 
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• 

K'g 

tidal-energy of the GC clue to its mass distribution 



inside the Galaxy gravity held 

• 

Vh 

potential energy clue to the Galaxy gravity held of the 



cluster barycenter 

• 

T c 

global kinetic energy of the GC 

• 

T' c 

intrinsic kinetic energy of the GC 

• 

J~<0 

kinetic energy of the cluster barycenter 

• 


self-potential energy of the GC 

• 

V' 

c 

Clausius’ virial energy of the cluster mass distribution 

• 

R 0 

mean orbit radius of the GC 

• 

M g (R 0 ) 

Galaxy mass fraction inside the sphere of radius R 0 

• 

act 

tidal radius at the maximum of Clausius’ virial 


in the two off-center spherical and homogeneous 
component model 
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